Vectorized Matlab codes for linear two-dimensional elasticity
نویسنده
چکیده
A vectorized Matlab implementation for the P1 finite element is provided for the twodimensional linear elasticity with mixed boundary conditions. Vectorization means that there is no loop over triangles. Numerical experiments show that our implementation is more efficient than the standard implementation with a loop over all triangles.
منابع مشابه
An efficient way to perform the assembly of finite element matrices in Matlab and Octave
We describe different optimization techniques to perform the assembly of finite element matrices in Matlab and Octave, from the standard approach to recent vectorized ones, without any low level language used. We finally obtain a simple and efficient vectorized algorithm able to compete in performance with dedicated software such as FreeFEM++. The principle of this assembly algorithm is general...
متن کاملTwo-Dimensional Elasticity Solution for Arbitrarily Supported Axially Functionally Graded Beams
First time, an analytical two-dimensional (2D) elasticity solution for arbitrarily supported axially functionally graded (FG) beam is developed. Linear gradation of the material property along the axis of the beam is considered. Using the strain displacement and constitutive relations, governing partial differential equations (PDEs) is obtained by employing Ressiner mixed var...
متن کاملVectorized Adaptive Quadrature in Matlab
Adaptive quadrature codes process a collection of subintervals one at a time. We show how to process them all simultaneously and so exploit vectorization and the use of fast built-in functions and array operations that are so important to efficient computation in Matlab. Using algebraic transformations we have made it just as easy for users to solve problems on infinite intervals and with moder...
متن کاملThe Nonlinear Bending Analysis for Circular Nano Plates Based on Modified Coupled Stress and Three- Dimensional Elasticity Theories
In this paper, the nonlinear bending analysis for annular circular nano plates is conducted based on the modified coupled stress and three-dimensional elasticity theories. For this purpose, the equilibrium equations, considering nonlinear strain terms, are calculated using the least energy potential method and solved by the numerical semi-analytical polynomial method. According to the previous ...
متن کاملON MAXWELL'S STRESS FUNCTIONS FOR SOLVING THREE DIMENSIONAL ELASTICITY PROBLEMS IN THE THEORY OF ELASTICITY
The governing equations of three dimensional elasticity problems include the six Beltrami-Michell stress compatibility equations, the three differential equations of equilibrium, and the six material constitutive relations; and these are usually solved subject to the boundary conditions. The system of fifteen differential equations is usually difficult to solve, and simplified methods are usual...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- Scientific Programming
دوره 15 شماره
صفحات -
تاریخ انتشار 2007